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Abstract: We compute power corrections to mean values of hadronic event shapes 
- the thrust and the C parameter — of tagged b quark events in electron positron 
annihilation, using the dispersive approach. We find that the leading power correc- 
tions are of the same type of 1/Q corrections as for event shapes in the massless 
case, with the same non-perturbative coefficient times a perturbatively calculable 
mass-dependent coefficient. The effect of the mass correction in the power correction 
is to reduce the latter by 10-30 % for tagged b events, for centre-of-mass energies 
ranging from the Z° peak down to 20 GeV. 
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1. Introduction 

In recent years the idea of an infrared-regular effective strong coupling at low scales 
(the strong coupling 'freezing' [1] or, more rigorously, the 'dispersive approach' of 
Dokshitzer, Marchesini and Webber [2]) has been employed for estimating hadroniza- 
tion corrections to various hadronic event shapes, using perturbative calculations 
[3-7]. These corrections arise in the form of power corrections of order 1/Q P , where 
Q is the centre-of-mass energy. Although the magnitude of the correction cannot 
be predicted by perturbative means, the power p and the relative coefficients among 
different observables can be calculated using the assumption of infrared freezing. 
This universality hypothesis has proved to give a phenomenologically fairly consis- 
tent picture of power corrections [8-13], and the same results could be derived using 
a different method, the so-called renormalon approach [14,15]. 

In this letter we explore a further check of the universality picture. Assuming 
the flavour independence of gluon radiation off quarks, the analytic structure of the 
strong coupling in the infrared, defined in an appropriate way, should not depend on 
the masses of the quarks, which radiate the gluon. Therefore, in the calculation of 
the power corrections to hadronic event shapes of tagged b events, the mass of the 
heavy quark enters in a completely perturbative manner. One may expect that the 
leading mass correction could be rrib/Q times the leading power correction (m& being 
the heavy-quark mass) with some calculable coefficient. If this coefficient is several 
times unity, then the mass corrections in the power corrections are expected to be 
significant enough to be a measurable effect. 
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A similar problem has already been considered by Nason and Webber in calcu- 
lating non-perturbative corrections to heavy-quark fragmentation in e + e~ annihila- 
tion [16]. Although the calculations are very similar, there is an essential difference 
in the results: for the fragmentation function there is a l/m b power correction, while 
for tagged events shapes there are mass corrections to the magnitude of the leading 
l/Q p power correction, which, however, vanish smoothly in the massless limit. 

Estimating the mass corrections in the power correction may be interesting from 
another point of view, too. Recently, next-to-leading order calculations of three- 
jet quantities in electron positron annihilation have been performed in which quark 
mass effects have been taken into account explicitly [17-19]. Using this theoretical 
input and the value of the 6-quark mass, the flavour independence of the strong 
interactions was demonstrated by determining the ratio of strong couplings, a b s /a^ ds 
[20,21]. One can turn around the argument and, assuming flavour independence, 
the 6-quark mass can be measured from the sensitivity of three-jet event shapes 
to mass effects [22-24]. Such a measurement requires an estimate of hadronization 
corrections. The traditional way of obtaining such estimates is by Monte Carlo event 
generators [25,26]. Those programs use the heavy-quark mass as input. As a result, 
the estimate of hadronization corrections brings a significant systematic error into 
the measurement of m&. For taking into account the hadronization correction in a 
different way, the simple formulae of power corrections presented in this paper can 
easily be incorporated in a fit to the 6-quark mass, which would hopefully result in 
a more stringent mass measurement. 

As mentioned above, the presence of the quark mass appears only in the per- 
turbative calculation. Therefore, the calculations that lead to the appearance of the 
Milan factor in the power corrections [5] can also be performed in the presence of 
the quark mass, resulting in a mass-dependent Milan factor. In this paper we make 
only the first step and calculate that part of power corrections to the mean value of 
two e + e~ event shapes, the thrust and the C parameter, which was termed 'naive 
contribution' in Ref. [5]. 

2. The dispersive approach 

In employing the dispersive approach of Ref. [2] one starts with assuming the validity 
of the following dispersion relation for the strong coupling: 

* s (k 2 ) = - J^-^p.tf) , p s {^) = -±Disc{a s (-^)} . (2.1) 

Further, it is also assumed that similarly to QED, the dominant effect of the running 
of a s on some QCD observable F may be represented in terms of the spectral function 
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p s (p 2 ) and a characteristic function J r (fJ 2 ): 

r°°dp 2 



F = a 8 (0) ^(0; ...)+/ Ps(p 2 ) ■ • •) (2-2) 



o p 



r^-Ps(p 2 ) [n^ 2 ; • • •) - ^(0; • • •)] , (2-3) 
Jo p A 1 



where we used Eq. (2.1) to express a s (0). The characteristic funtion JF(/i 2 ; . . .) is a 
function of dimensionless ratios of e = p 2 /Q 2 , p = 4m 2 /Q 2 , where m b is the heavy- 
quark mass and {y}, the collection of any further relevant dimensionless parameters, 
e.g. the jet shape S in our present considerations. We obtain T by computing the 
one-loop graphs, corresponding to the physical process under consideration, with a 
non-zero gluon mass e ^ 0, and dividing by a s : 

T{e- p,S)= J d$({^}; p, e)M({xi}; p, e) S({x t }- p, e) , (2.4) 

where d$({xj};p, e) denotes the phase space in terms of the independent variables 
M. is the proper squared amplitude divided by a s , and S stands for the event 
shape variable. Introducing the effective coupling a e s(p 2 ), defined in terms of the 
spectral function by 

PV)=/^, (2.5) 



we can integrate Eq. (2.3) by parts to obtain 

F = r ^ a cS (p 2 ) T(p 2 ) = [°°- a cS (Q 2 e) P{e) , (2.6) 
Jo a z Jo e 



where 

T{e) = -e^ . (2.7) 

Using the definition Eq. (2.5) and the dispersion relation Eq. (2.1), we can deduce 
that in the perturbative domain, q s < 1, the standard and effective couplings are 
approximately the same [2]: 

a cS (p 2 ) = a s (p 2 ) + 0{a 3 s ) . (2.8) 

Thus we may interpret a c fj as an effective coupling that extends the physical per- 
turbative coupling into the non-perturbative domain. If the effective coupling has a 
non-perturbative component 5a e s(p 2 ), with support limited to low values of p 2 , the 
corresponding contribution to F, 

5F = J — 5a eS (Q 2 e)F(e), (2.9) 

will have a Q 2 dependence determined by the small-e behaviour of T . 
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We shall be interested in the leading power behaviour of 5F, which will be called 
the leading power correction to F, 

5F = F pow + 0(l/Q p ) . (2.10) 

In order to determine F pow , we recall that power-suppressed contributions to 5F can 
arise from only those terms in T that are non-analytic at e = [2]. Therefore, the 
leading power correction _pp° w will be obtained from the leading non-analytic term 
at e = in JF, which we denote by Cf /2-k. 



3. Calculations 



In calculating the function bT for the e + e~^QQg event shapes, we find three sources 
of quark-mass dependence: (i) restriction of phase space; (ii) mass corrections in the 
matrix element; (iii) mass corrections in the definition of the event shape. The double 
differential cross section for the production of a massive quark antiquark pair and a 
massive gluon given in terms of the scaled quark, antiquark and gluon energies, x, x 
and x g = 2 — x — x was derived in Ref. [16]. This cross section is an analytic function 
of the gluon mass and, therefore, in calculating SJ-, the e-dependent terms can be 
dropped completely. In the case of the vector current contribution the cross section 
for zero gluon mass is given by 



1 do\ 



oy da; dx (3 2n 



-\{x 2 



p, x 2 



8 (1 - x) 2 (1 - x) 2 
1 / 1 



(3.1) 



+ 



X 



x) 2 



for the axial vector current contribution we have 



1 da, 



a a dx dx (3 2ir 



+ 



-\{x 2 -p, x 2 -p, x 2 g ) 
8 (1 - x) 2 {1 - x) 2 
1 ( 1 



(3.2) 



+ 



x) 



x? + 2 



X) 



X) 



where (3 = a/1 — p is the quark velocitey and X(x, y, z) 
2zx. In Eqs. (3.1) and (3.2) 



x 2 + y 2 + z 2 — 2xy — 2y z ■ 



a v = a (l + \ p) (3 , 



(3.3) 



and 



a A = a (l-p) (3 (3.4) 

are the Born cross sections for heavy-quark production by a vector and an axial 
vector current respectively, o"o being the massless quark Born cross section. The 
corresponding phase space was also given in Ref. [16]: 

d$(x, x;e,p) = dx dx 0^ — | X(x 2 — p,x 2 — p, x 2 — Aefj (3.5) 
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The phase-space boundary in this equation gives x_ < x < x + , where 



(2- x )(l- x -e + ip)± s /(T'-p)[(l-x-er-ep] 

X± = 2(l-s) + i„ (3 ' 6) 

and 

y/p < x < 1 — e — ^fep . (3.7) 

For large values of the gluon momentum the phase-space boundary does not yield 
leading non-analytic contributions in e. Therefore, the leading non-analytic contri- 
bution to comes from the soft-gluon emission region [2,16], where both the matrix 
elements, in Eqs. (3.1) and (3.2), and the phase space, in Eq. (3.5), can be expanded 
in x g = y + y ~ \fe. This expansion for the cross section formulae results in 



1 day a Cf a 



a VtA dydy (3 2% 



2 p(y + y) 2 



yy 2 y 2 y 2 



(3.8) 



both for vector and axial vector current contributions, where y = 1 — x and y = 1—x. 
The lower boundary of the phase space in x g = y + y ~ -Je region is approximated 
as follows: 

-\ K* 2 - P^ 2 - P,x 2 g ~ 4e) » 2yy - P - (y + yf - 2e (1 - p) = . (3.9) 

Thus, to obtain bT for mean values of event shapes, we have to calculate the following 
integral 



ST £ d 



J dy dyQ(y - Jep) Q(y - Jep) 6(1 -p/2-y-y) (3.10) 
x&(4yy-p(y + y) 2 -4e(l-p)) 



2 p (y + y) 2 \ _ 

x— 1-t — )S{y,y;p,£ ). 

yy \ 4 yy J 

In this equation, the exact form of the upper boundary does not influence those 
non-analytic terms that give the leading power correction. 

The mass corrections in the physical quantity introduce an observable-dependent, 
but perturbatively calculable mass dependence. As simple application of formula 
(3.10), we consider two event shapes, the C parameter and the thrust T. The C 
parameter [27] is derived from the eigenvalues of the infrared-safe momentum tensor 

^' = £#/£bU (3.ii) 

a \Pa\ ' a 

where the sum on a runs over all final-state hadrons and p z a is the ith component of 
the three-momentum p a of hadron a in the cm. system. The tensor 9 is normalized 
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to have unit trace. In terms of the eigenvalues Aj of the 3x3 matrix 6, the global 
shape parameter C is defined as 



C = 3 (AiA 2 + A 2 A 3 + A 3 Ai) . (3.12) 

The thrust [28] is given by 

rp ( J2a \Pa ■ n T \ /Q 1 Q\ 

1 = max r-^-j — , (3.13) 



where the sum runs over all final-state particles and the thrust axis Ht is chosen 
to maximize the expression. For three partons in the final state, the thrust can be 
written as [29] 

rp = 2 max{|pi|, \p 2 \, \p 3 \} ^ 

\Pa\ 

Instead of the thrust, we shall consider its deviation from unity, t = 1 — T. 

For these two event shapes the contribution of the non-perturbative gluon will 
just add to that of the underlying perturbative event (contribution of many soft per- 
turbative gluons), and the power correction based on the presence of just a single 
non-perturbative gluon will remain valid [13]. This argument can be used indepen- 
dently of the mass of the leading quark pair. As discussed in Ref. [5], the 'naive 
contribution' to the mean value of the C parameter is obtained from Eq. (3.12) 
in the soft-gluon approximation, with gluon mass set to zero. The corresponding 
formula with non-zero quark masses is the following: 

(3 s y + y \ 4 yy J 

Similarly, we need only the leading term in the soft-gluon expansion of the function 
t(y, y; p) with the gluon mass neglected, 

t(y, mp)~ \ (4== - t~ — \(v -y)<3>(y-y) + (y-y)<3>(y-y)}) ■ (3.16) 
2 vv 1 - p 1 - p L V 

To perform the integrations in Eq. (3.10) we introduce polar coordinates, y = 
r cos 0, y — r sin 0. In terms of the variables r and <fi 

d$(e) = r dr d0 0(0 - 5) 0(|tt - 5 - 0) 0(r+(0) - r) 0(r - r_(0, e)) , (3.17) 

where 



(l + sin2^)- 1 2 g (l-sin2£) p 

r + (0) — . . ==r . r_ 0,5) = W — -. sin 2d = . 3.18 

+ VI + sin 20 ' ; V sin 20 - sin 25 ' 2-p v ; 
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The r integral is trivial and the integral over <p can be performed after making the 
shift — > — 7r/4. After performing the differentiation in Eq. (3.10), for the C 
parameter, S = C, we obtain 

= 6irVi(l-|v^ + 0(p)) . (3.20) 

From the expansion at small values of p, we see that our formula reproduces the 
known zero mass result [1]. In the case of thrust, S = t, one finds 



V p V i -p v 



(3.21) 



where 

2 -p± \f2 . 
a± = arctan — , 3.22) 

JpO^J) 

Expanding in p we find 

8F® = Ay/i (l-yv^+Olp)) , (3-23) 

which shows that the apparent \Je/ p behaviour in Eq. (3.21) is in fact a y/e be- 
haviour with multiplicative mass corrections that are regular for vanishing quark 
mass. Setting p = 0, we find agreement with the known zero-mass result [1]. 

In Fig. 1 we plotted the functions jF(e) in Eq. (2.7), obtained from numerical 
integration for representative values of p, with T as given in Eq. (2.4), but with S 
taken in the soft-gluon approximation (Eq. (3.15) for S = C in Fig. la and Eq. (3.16) 
for S — t in Fig. lb). In the same figure, the solid lines show the analytic results 
5J-( C \ Eq. (3.19), and 8J-^\ Eq. (3.21). We clearly see the leading ^fe behaviour for 
small values of e together with the correct quark-mass dependence. 

Having calculated the derivative of the characteristic functions for the various 
event shapes, we can use Eq. (2.10) to obtain the power correction F pow . Following 
Ref. [2], we introduce the moment integral 

M P = ^ -e p 8a cS (eQ 2 ). (3.24) 
Ztt Jo e 

In terms of this non-perturbative parameter, the 'naive contribution' to the power 
corrections to the mean value of C and t are given by the p = \ moment of 8a c s as 



{ C )v™ = 6iv^- 1-3^ + ^ (3.25) 
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Figure 1: Derivative of the characteristic function for the mean value of the (a) C param- 
eter (b) 1 — T for values of p = 0, 0.09, 0.36, 0.81. 



and 

with the same value of the moment integral as phenomenologically deduced from 
untagged samples. To obtain these leading mass corrections, valid for small values of 
the tagged-quark mass, we used Eqs. (3.20) and (3.23). These results show that at 
LEPI energies the mass correction reduces the magnitude of the power correction by 
about 10% for the mean value of the C parameter and by about 7% for the mean 
value of t (with m b taken to be the MS mass at the centre-of-mass energy). 

A two-loop analysis of power corrections to event shapes in the massless-quark 
case revealed that there is some freedom in the definition of the 'naive contribution'. 
With a different definition for this term, the other two ('inclusive' and 'non-inclusive') 
contributions change in such a way that the sum of the three terms remains unam- 
biguous. The universal result can be summarized by a simple multiplicative correc- 
tion factor, the so-called Milan factor [5]. 1 A two-loop calculation with mass effects 
taken into account would modify our results in Eqs. (3.25) and (3.26) in two ways. 
On the one hand, these equations would acquire an overall Milan factor of 1.5 and on 
the other the coefficients in front of the m b /Q corrections would change. For the rest 
of this paper, we shall neglect the latter modification, but shall include the Milan 
factor. 



4. Merging perturbative and non-perturbative contributions 

Once we have the leading power correction, we have to combine it with the per- 
turbative prediction for the same physical quantity to obtain the full theoretical 

1 Recently , it was pointed out that the correct value of the Milan factor is A4 « 1.5 [30]. 
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prediction: 

F{Q) = F pert (Q) + F pow (Q). (4.1) 

The two contributions are separately ill-defined because the perturbative part is 
given by an expansion that is factorially divergent, while in the power correction, 
the non-perturbative part of the effective coupling Sa e g contains an ill-defined all- 
order subtraction of the pure perturbative part off the full effective coupling a e s. 
The sum of the two contributions is finite. At fixed order in perturbation theory 
they should be merged in such a way that the terms that would grow factorially 
in an all-order result cancel order by order. We follow the prescription of Ref. [5] 
for this merging, where the moment integral A\ was approximated with the same 
integral up to an infrared scale /xj, 

A 1 ~ ^% r'dk \a s (k) - a p s eit (k)] . (4.2) 

71 JO L J 

In this equation the integrand is the non-perturbative component of the strong cou- 
pling, i.e. the difference of the strong coupling a s as given by Eq. (2.1) and the 
perturbative coupling a Pert . Above the infrared scale fi\, this non-perturbative cou- 
pling is assumed to give negligible contribution to the moment integral (see Eq. (2.8) 
and the following paragraph). 

The integral of the strong coupling in Eq. (4.2) can be expressed in terms of a 
phenomenological parameter a , 

/•« 

dka s (k) = /iia (iii) . (4.3) 

/o 

The value of this parameter depends on the infrared matching scale. For pii — 2 GeV, 
fits of event shapes of untagged events gave a ~ 0.5 [1,9,11-13]. 

To calculate the integral of the perturbative coupling a pcrt , we use its one- loop 
expression given as a geometric series: 

where (3^ is the one-loop beta function 

^ Nt) = ^C A -\N u (4.5) 

for Nf — 4 light fermion flavours, and a pert (<5) is a four-flavour perturbative coupling 
in the physical (CMW) renormalization scheme [31]. 

For the practically interesting cases the perturbative prediction for the event 
shape the perturbative prediction _F pert is known to second order in a s , 

F p ert(g . 2) = ^(Q) B(m b , Q) + (^^] C(m b , Q) + O ((a? 5 (Q)) 3 ) , (4.6) 



i: 
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with Born and correction coefficients, B(rrib, Q) and C(rrib, Q) respectively, calculated 
in the MS scheme at scale Q, and a^ s (Q) is the four-flavour strong coupling defined 
in the MS renormalization scheme. In such cases the summation in Eq. (4.4) should 
be truncated at O(a^), that is at £ = 1. Then the integral of the perturbative 
coupling in Eq. (4.2) gives 




«r rt (Q) 



(4.7) 



In order to use the MS coupling everywhere in the final result, one has to make the 
shift _ 

< rt -« s MS (l + ^fH (4-8) 
in Eq. (4.7), with K defined as 

* = Cf * (is Nt = i - (49) 

The coupling in the definition of the non-perturbative parameter a® is the physical 
coupling; therefore, the value of a® does not depend on the chosen scheme. 

Collecting the various contributions, for the event shape S(— t, or C) we find 

F {S) (Q) = a s B + a 2 s C (4.10) 



7r y/e Q \ 



a, — a 2 




+ K 



where a s = a^ LS (Q)/2n for Nf — 4 flavours and «o = ao/2n. In order to use the 
standard Nf = 5 flavour MS strong coupling, we have to make the shift 



a s -> a s + -a 2 s In — f . (4.11) 




We could make a more sophisticated description of the heavy-quark threshold in the 
running coupling (see Refs. [32,33]), but the usual centre-of-mass energy for event 
shapes being far from the threshold, the difference is negligible. For instance, using 
the prescription of Ref. [33] changes the physical prediction F^ S \Q) by less than 
1 % at Q — 20 GeV and by about 0.1 % at the Z° peak. The renormalization scale 
dependence of this result can be studied in the usual way. To show the effect of the 
mass correction, we plot for S = C in Fig. 2a and for S — t in Fig. 2b for 
/ii — 2 GeV, a = 0.5 and the world average of the strong coupling at the Z° peak, 
a s (M z ) = 0.119 [34]. The solid lines show the results for the central value of the MS 
6-quark mass at the given hard scattering scale run from m fe (m b ) = 4.3 GeV [32]. The 
dotted line is the next-to-leading order perturbative prediction, with mass effects 
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included, and the dashed lines represent the result when mass effects are present 
in the perturbative prediction plus power corrections without mass effects. The 
perturbative coefficients B and C were obtained using the zbb4 program [35]. We 
can observe that the mass effect in the power correction in tagged b samples is 
important for centre-of-mass energies below about 45 GeV. 



0.15 



F""+F p ""(m b =0) 

F P« +F PC»: (mb=4 3 Gf , v) 



ppe« 

F p«, 1+F p«„ (mb=0) 

F P«t +F l»w (mt)= 4 3 GeV) 



a 



0.06 — 



0.05 



10 20 30 40 50 60 70 



90 100 110 



0.04 



10 20 30 40 50 60 70 

Q 



90 100 110 



Figure 2: Mean value of (a) the C parameter and (b) 1 — T, where T is the thrust. Dotted: 
next-to-leading order perturbative result, dashed: perturbative + power correction without 
mass effect, solid: perturbative + power correction with mass effect. 



If we compare our results for the power corrections of tagged b events to the 
hadronization corrections obtained using Monte Carlo programs, we find significantly 
smaller corrections from our model. For instance, the double ratio H$ = / h\ S \ 
where is the hadron level value divided by the parton level value for event 

shape S and for the quark flavour x (heavy, or light), was determined in Ref. [24] 
for events in Z° decays using JETSET parton shower model together with the Lund 
string fragmentation model. In this model the mass effects are introduced only by 
kinematic constraints to the phase space at each parton branching in the shower 
evolution. Table 1 shows the values obtained for H^J using the two models. 

The values for the double ratios obtained 
from the dispersive approach are very close to 
1, indicating that the relative power corrections 
are very similar in the 6-tagged and w<is-tagged 
samples. 2 The situation is very different for the 
JETSET estimates. The results, when trans- 
lated for single ratios of hadronization correc- 
tions on the tagged sample Hj; S \ mean that 
the hadronization correction from the disper- 
sive approach is about half of that from JET- 



s 


JETSET 


Disp. model 


c 
t 


1.175 
1.142 


1.000 
1.001 



Table 1: Comparison of hadroniza- 
tion corrections obtained using JET- 
SET and from the dispersive model 
via the double ratios and H^J V 



2 This result depends very weakly on the value of the Milan factor. 
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SET. It would be interesting to learn whether or not this difference changes if the 
mass effects are taken into account in the dynamics of the showering in JETSET, or 
it is rather due to the different secondary decays of heavy hadrons. The latter case 
may undermine the usefulness of estimating hadronization corrections from power 
corrections for tagged b events. 

5. Conclusion 

In this paper we have calculated the 'naive contribution', as defined in Ref. [5], for the 
power corrections of 6-tagged event shapes thrust and C parameter in e + e~ annihila- 
tion. At LEPI energy the explicit analytic formulae predict about 7-10 % reduction 
of the hadronization correction for the shapes obtained from 5-tagged samples with 
respect to the untagged case. This reduction effect will be even more important for 
t-tagged samples at the NLC. For instance, for m t = 175 GeV and Q = 500 GeV, 
we expect the power correction to be about two thirds of the corresponding massless 
case for the C parameter and t. 

We also presented predictions for the mean values of the event shapes by merging 
the perturbative and non-perturbative contributions. The only non-perturbative 
input in the prediction is the non-perturbative parameter a , taken from fits of 
theoretical predictions for the mean values of event shapes to those of light primary 
quark samples. We found that the mass effect in the power correction for tagged 
b samples is not too profound for centre-of-mass energies above about 45 GeV and 
is significantly smaller than the hadronization correction predicted by the current 
version of JETSET, which does not take into account the heavy-quark mass in the 
dynamics of fragmentation, but includes decays of heavy hadrons. The results may 
be used directly in heavy-quark mass measurements. 

In our view the flavour independence of the strong coupling is a sufficiently weak 
assumption so that if it is feasible with the dispersive approach to calculate non- 
perturbative effects in hadronic event shapes, then the primary quark mass effects 
should be perturbatively calculable. Therefore, our results, once confronted with 
experiment, can serve as a further check of perturbatively calculable hadronization 
corrections. 
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